AD-A196 


UWLLAbbiHtU 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whnn  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


report  number 


2.  GOVT  accession  no 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CAT  ALOG  NUMBER 


kjlV|»av»V|f|V*|V 


AFIT/CI/NR  88-  gS 


4.  TITLE  (and  Subtitle)  ^lJHo0  fi  ^  CjTAuJ/OG  £MplftlCAL 

ColWl-fCuArioOS  POA  PG.£OicT;PG  C.W.AC*  PlLt>pf\GfiT(oiJ 
f1°  A  i3u(UkJ  IkJG  5<3UD  PCoftuAWT  CRAifU 


5.  type  of  report  a  PERIOD  covered 


MS  THESIS 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR(a) 


0  £F.p  R.Ly  A.  NjMi3 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

AFIT  STUDENT  AT:  PtkJM  JyL  i/AwjA  5TATL 

upi  vfi.K.5  try 


CONTROLLING  OFFICE  NAME  AND  ADDRESS 


14.  MONITORING  AGENCY  NAME  &  ADDRESS  (If  different  from  Controlling  Office; 

AFIT/ NR 

Wright-Patterson  AFB  OH  45433-6583 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


8.  CONTRACT  OR  GRANT  NUMBERftJ 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

1988 


13.  NUMBER  OF  PAGES 

S9 


15.  SECURITY  CLASS,  (of  this  report) 


UNCLASSIFIED 


1  Sa.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


DISTRIBUTED  UNLIMITED:  APPROVED  FOR  PUBLIC  RELEASE 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  It  different  from  Report)  * 


SAME  AS  REPORT 


10.  SUPPLEMENTARY  NOTES 


DT!C 


AUG  0  4  1988 


Approved  for  Publ ic^Release: \  IAW  AFR  190-1  < 

LYNN  E.  WOLAVER  iyJJ^vs  “'fA  W 

Dean  for  Research  gxd  Professional  Development  v 
Air  Force  Institute  of  Technology 


i?iaEinagwngpnwja«ti;KiMiicicga.f:fci 


19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 


20.  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

ATTACHED 


1  JAN  73  1473  EDITION  OF  I  NOV  65  IS  OBSOLETE 


uricLAssiFef 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wien  Data  Entarad) 


m 


WA'^W'i i 


ABSTRACT 


The  importance  of  crack  propagation  in  solid  rocket  motors  is 
widely  recognized.  However,  the  processes  of  crack  propagation  and 
branching  in  burning  solid  propellants  are  not  as  yet  well  understood. 
These  processes  could  be  instrumental  in  creating  large  specific 
surface  areas  for  burning  causing  rocket  motor  failure.-^ Numerous 
samples  have  been  recovered  from  interrupted-burning  experiments. 
Depending  on  the  pressurization  rate  and  sample  geometry,  different 
modes  of  crack  propagation  and/or  branching  were  observed.  This 
problem  was  analyzed  through  the  use  of  basic  physical  principles, 
current  theories  on  crack  propagation  in  viscoelatic  media,  and 
experimental  results.  \^A  set  of  governing  dimensionless  parameters 
which  control  and  characterize  the  degree  of  damage  was  obtained.  ^-Key 
factors  influencing  the  crack  combustion,  propagation  and  branching 
process  considered  in  the  dimensional  analysis  include  chamber 
pressurization  rate,  initial  crack  length,  initial  temperature,  sample 
geometry,  and  mechanical  and  thermal  properties  of  the  solid 

i 

propellant.  .  Based  on  experimental  results,  a  multiple  regression 
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analysis  was  performed  to  obtain  a  set  of  correlations  between  the 
dimensionless  parameters.  The  procedures  developed  in  this  work 
demonstrate  a  methodology  for  obtaining  similar  correlations  which 
could  be  used  as  guidelines  for  safe  operating  conditions  of  rocket 
motors. 
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ABSTRACT 


The  importance  of  crack  propagation  in  solid  rocket  motors  is 
widely  recognized.  However,  the  processes  of  crack  propagation  and 
branching  in  burning  solid  propellants  are  not  as  yet  well  understood. 
These  processes  could  be  instrumental  in  creating  large  specific 
surface  areas  for  burning  causing  rocket  motor  failure.  Numerous 
samples  have  been  recovered  from  interrupted-burning  experiments. 
Depending  on  the  pressurization  rate  and  sample  geometry,  different 
modes  of  crack  propagation  and/or  branching  were  observed.  This 
problem  was  analyzed  through  the  use  of  basic  physical  principles, 
current  theories  on  crack  propagation  in  viscoelatic  media,  and 
experimental  results.  A  set  of  governing  dimensionless  parameters 
which  control  and  characterize  the  degree  of  damage  was  obtained.  Key 
factors  influencing  the  crack  combustion,  propagation  and  branching 
process  considered  in  the  dimensional  analysis  include  chamber 
pressurization  rate,  initial  crack  length,  initial  temperature,  sample 
geometry,  and  mechanical  and  thermal  properties  of  the  solid 
propellant.  Based  on  experimental  results,  a  multiple  regression 
analysis  was  performed  to  obtain  a  set  of  correlations  between  the 
dimensionless  parameters.  The  procedures  developed  in  this  work 
demonstrate  a  methodology  for  obtaining  similar  correlations  which 
could  be  used  as  guidelines  for  safe  operating  conditions  of  rocket 
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Chapter  1 


INTRODUCTION 


1 . 1  Motivation  and  Objectives 

Solid  propellant  rocket  safety  is  a  key  factor  in  several  of 

today's  most  important  civilian  and  military  propulsion  systems.  With 

both  civilian  and  military  agencies  presently  pursuing  solid  rocket 

booster  safety  and  ballistic  risk  assessment  programs,  near-term 

advances  in  this  area  are  foreseeable.  Of  particular  concern  is 

overpressurization  and  subsequent  catastrophic  failure  of  a  rocket  motor 

due  to  anomalous  burning,  which  may  be  caused  by  the  presence  of  voids, 

cracks,  or  debonded  regions  in  the  propellant  grain.  Such  regions  are 

susceptible  to  further  damage  under  transient  loading  conditions.  A 

theoretical  model  has  been  developed  to  predict  the  flame  spreading  and 

combustion  processes  inside  a  propellant  cavity1  and  various  models  have 

also  been  developed  for  crack  initiation  and  growth  in  viscoelastic, 
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inert  material.  However,  the  coupling  effects  of  these  two  processes 
are  complicated  and  not  totally  understood.  Based  on  previous  work  by 

6_g 

Kuo  and  co-workers  on  crack  growth  in  burning  propellant  grains,  a 
set  of  governing  dimensionless  parameters  have  been  identified  to  obtain 
empirical  correlations. 

The  primary  objectives  of  this  research  are: 

a)  to  define  a  group  of  meaningful  dimensionless  parameters  for 
characterizing  the  operating  conditions  and  resulting  damage  of  a  solid 
propellant  sample; 

b)  to  obtain  test  data  necessary  for  developing  correlations 
between  the  dimensionless  parameters;  and, 


"\v; 


c)  to  establish  correlations,  demonstrating  a  methodology  for 
generating  useful  safety  criteria. 


1.2  Brief  Survey  of  Related  Work 


Early  studies  of  the  crack  propagation  phenomena  under  burning 
conditions  by  Siefert  and  Kuou  provided  good  insight  into  the  basic 
mechanisms  involved.  They  studied  the  crack  propagation  process  of  a 
burning  propellant  sample  of  constant  geometry  under  varying 
pressurization  rates.  Experimental  evidence  supported  the  assumption 
that  the  crack  propagation  velocity  was  limited  by  the  rate  of 
deformation  of  the  crack  cavity  walls.  In  this  case,  a  critical  stress 
state  was  maintained  in  the  failure  zone  at  the  crack  tip.  A  nearly 
linear  relationship  between  the  crack  velocity  and  pressurization  rate 


was  observed, 


V  (m/s)  =  0. 086 [dp/dt( MPa/s)] 


0.637 


For  low  pressurization  rates,  the  crack  tip  maintained  its  original 

contour.  However,  at  pressurization  rates  on  the  order  of  20  GPa/s,  the 

crack  tip  appeared  to  square  off  and  a  fan  structure  developed 

suggesting  the  occurrence  of  crack  branching  and  microstructure  damage. 

It  is  interesting  to  note  that  the  maximum  crack  velocities  measured 

correspond  to  terminal  velocities  of  inert  propellants  and  rubber 

9 

reported  by  Gent  and  Marteny. 


Further  studies  by  Kuo,  Moreci  and  Mantzaras  ’  have  identified 
four  different  modes  of  damage  which  were  related  to  the  pressurization 
rate  of  the  crack  cavity.  Figure  1  shows  the  measured  crack  propagation 
velocity  as  a  function  of  pressurization  rate  and  where  each  of  these 
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modes  was  observed.  Mode  A  represents  little  or  no  increase  in  crack 
length  and  corresponds  to  very  low  pressurization  rates.  No  departure 
from  normal  burning  is  expected  for  this  case.  Mode  B  represents  the 
propagation  of  a  single  crack.  The  rate  of  crack  growth  increases 
almost  linearly  with  pressurization  rate  in  the  range  of  1.4  to  15.0 
GPa/s.  Although  burning  is  enhanced  at  the  crack  tip,  the  overall 
increase  in  burning  due  to  the  damage  is  not  significant.  Mode  C 
represents  an  initial  single  crack  propagation  with  local  branching  at 
various  locations.  This  was  observed  for  intermediate  pressurization 
rates  and  has  the  potential  of  significantly  increasing  the  burning 
surface  area.  Mode  D  represents  crack  branching  from  the  initial  crack 
tip  location  in  either  discrete  macrocracks  or  a  fan  of  numerous  smaller 
cracks.  This  type  of  damage  was  characteristic  of  higher  pressurization 
rates  and  also  can  significantly  increase  the  burning  surface  area. 

Some  of  the  data  from  these  experiments,  as  well  as  newly  acquired  data, 
were  used  in  the  current  analysis. 
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Chapter  2 


EXPERIMENTAL  SETUP 


2. 1  Selection  of  Sample  Propellant  Type 


In  selecting  a  propellant  sample  for  this  research,  the  following 
criteria  were  considered: 

a)  the  propellant  must  be  available  for  a  basic  study  in  a 
university  environment; 

b)  its  combustion  behavior  must  be  well  characterized; 

c)  its  mechanical  properties  and  failure  criterion  under  non-buring 
conditions  must  be  known;  and, 

d)  propellants  of  the  same  family  must  be  available  so  that  similar 
tests  can  be  conducted  to  determine  the  effect  of  burning  rate  and 
oxidizer  size. 

Two  composite  propellants  (ARC  5051  and  ARC  4525)  were  selected. 
Each  of  these  contain  27%  HTPB  and  73?  AP  with  average  oxidizer  particle 
sizes  of  200um  (ARC  5051)  and  20ym  (ARC  4525).  Due  to  the  similarities 
in  composition  thermodynamic  properties  such  as  enthalpy,  density, 
burning  rate,  etc.,  are  nearly  equal.  Their  combustion  behavior  and 
mechanical  properties  have  been  well  characterized.  In  addition, 
propellants  of  similar  compostion  but  with  different  burning  rates  are 
available  for  future  studies.  Although  these  propellants  have  a  lower 
solids  loading  than  the  typical  high-energy  solid  propellants  used  in 
most  of  today's  advanced  rocket  propulsion  systems,  once  the  basic 
method  for  determining  crack  combustion  correlations  has  been 
established,  similar  studies  can  be  made  for  other  propellants. 
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2.2  Sample  Configuration 

The  two-dimensional  sample  configuration  used  in  previous  studies 
is  adopted  for  this  study  (Fig.  2).  This  general  geometry  is  similar  to 
one  point  of  a  star  grain  or  to  a  segment  of  a  grain  containing  a  void. 
This  configuration  provides  an  initial  crack  of  variable  length  and 
propagation  along  the  horizontal  axis.  The  geometric  parameters  i|i  and 
tp.  can  also  be  varied  to  represent  different  local  boundary  conitions  of 
the  damaged  zone. 


2. 3  Test  Setup 

The  test  setup  was  designed  to  induce  combustion  and  mechanical 
damage  in  a  localized  zone  of  a  pre-cut  propellant  sample  by  subjecting 
it  to  high-pressure,  high-temperature  gases  and  to  interupt  the  buring 
by  rapid  depressurization  for  sample  recovery.  In  these  tests 
variations  of  several  independent  parameters  were  made,  including 
pressurization  rate,  initial  crack  length,  sample  geometry  and 
propellant  composition.  The  damage  was  then  measured  in  terms  of  crack 


propagation  velocity,  extended  crack  length  and  number  of  cracks. 


2.3.1  Test  Chamber 

The  test  rig  (Fig  3)  consisted  of  a  high-pressure,  windowed  chamber 
with  the  following  subsystems: 

a)  a  driving  motor  to  pressurize  the  free  volume  of  the  chamber; 

b)  a  rupture  disk  and  venting  port  to  rapidly  depressurize  the 
chamber;  and, 

c)  a  nitrogen  injection  system  to  cool  the  exposed  surfaces  and 
reduce  the  chance  of  sample  reignition. 


Test  Chamber  Schematic 
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2.3. 1.1  Driving  Motor 

The  driving  motor  contained  four  key  components  (Fig  4).  The 
ignition  primer  (type  M52A3B1),  when  supplied  with  three  to  four  volts 
DC,  would  fire  and  in  turn  ignite  the  propellant  shavings.  The  shavings 
burn  extremely  fast  due  to  their  very  high  burning  surface  area.  Their 
purpose  was  to  supply  hot  gases  ad  boost  the  pressure  in  the  segment 
containing  the  main  igniter  charge.  The  igniter  charge  would  then  burn, 
exhausting  through  the  perforated  nozzle  into  the  main  chamber.  By 
varying  the  mass  of  the  main  igniter  charge,  pressurization  rates 
ranging  from  near  zero  to  50  Gpa/s  were  achieved. 


2.3. 1.2  Rupture  Disk  and  Venting  Port 


Igniter 

Charge 


Shavings 


injected  following  depressurization  to  cool  the  surfaces  of  the  free 
volume. 

2.3.2  Data  Acquisition  System 

The  data  needed  from  these  tests  to  calculate  key  parameters 
consisted  of  a  pressure-time  trace  of  the  crack  cavity  as  well  as  a 
visual  record  of  the  crack  itself  as  it  propagated.  A  complete 
schematic  of  the  data  acquisition  system  is  shown  in  Fig  5. 

2. 3. 2.1  Pressure  Measurements 

The  pressure  measuring  device  selected  was  a  miniature 
piezoelectric  quartz  transducer  (Kistler  Model  60 1 B 1 ) .  This  transducer 
has  a  10%  to  90%  rise  time  of  3.0  microseconds  over  the  range  of  0  to 
15,000  psi.  This  was  more  than  adequate  for  measuring  pressurization 
rates  of  up  to  50  GPa/s  with  rise  times  on  the  order  of  1  millisecond. 

To  protect  the  transducers  from  the  high  temperature  product  gases,  they 
were  recess-mounted  in  water-cooled  adapters.  An  insulated 
high-impedence  cable  with  a  nominal  capacitance  of  30  pF/ft  was  used  to 
transmit  the  signal  from  the  transducer  to  the  charge  amplifiers.  The 
signal  attenuation  under  the  test  conditions  was  less  than  one  percent. 

The  charge  amplifiers  (Kistler  Model  504E)  output  the  transducer 
charge  signal  as  an  amplified  voltage.  Calibration  is  achieved  in  the 
charge  amplifier  by  supplying  a  DC  voltage  of  equal  magnitude  to  that  of 
the  pressure  transducer  sensitivity.  The  output  voltage  of  the  charge 
amplifier  is  then  adjusted  to  satisfy  the  equation 


-V.  x  1 OOOpF 

y  _  - LD - 

out  ”  Range  x  Sensitivity 


(2) 


A  typical  range  setting  of  2000  psi/volt  resulted  in  a  calibration 
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DIGITAL 
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output  of  0.5  volts.  A  long  time  constant  was  used  to  provide  a 
continuous  signal. 

The  output  from  the  charge  amplifiers  was  recorded  on  a  Physical 
Data  Model  51 5A  Transient  Waveform  Recorder.  It  recorded  the  voltage 
signals  at  sampling  intervals  of  5  to  50  microseconds  for  the  initial 
pressurization  period  and  then  shifted  to  a  sampling  rate  of  2.5  to  5 
milliseconds  in  order  to  extend  the  total  time  recorded.  The  limiting 
factor  was  its  maximum  capacity  of  4096  data  points.  The  recording 
mechanism  was  triggered  by  a  small  signal  from  the  transducer  in  the 
crack  cavity.  A  negative  delay  of  512  data  points  was  used  to  insure 
the  beginning  of  the  event  was  recorded.  The  recorded  signal  could  then 
be  played  back  to  a  Nicolet  Digital  Oscilloscope  for  analysis  and 
recording  on  floppy  disk  or  to  an  X-Y  plotte”  for  a  hardcopy. 

2.3. 2. 2  High  Speed  Film  Recording 

A  Hycam,  16  mm  camera  was  used  to  record  the  ignition,  combustion, 
and  crack  propagation  of  the  propellant  sample  through  a  veiwing  port  in 
the  chamber.  Typical  framing  rates  were  5,000  frames  per  second  in  a 
half  frame  format  resulting  in  10,000  pictures  per  second.  In  order  to 
insure  the  event  occurred  before  the  roll  of  film  was  completely 
exposed,  the  camera  was  used  to  trigger  the  ignition  of  the  driving 
motor.  The  camera  was  started  and  when  a  prescribed  footage  of  film  was 
exposed,  allowing  the  camera  to  reach  the  set  film  speed,  an  internal 
switch  would  close  the  igniter  circuit.  This  would  also  trigger  a 
common  time  pulse  which  was  recorded  both  on  the  film  and  on  the 
Physical  Data  Recorder.  Another  light  pulse  with  a  frequency  of  1  KHz 
was  also  recorded  on  the  film,  to  one  side  of  the  main  image,  in  order 
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to  verify  the  framing  rate. 

2.4  Test  Procedure 

Appendix  A  contains  a  complete  checklist  and  sample  data  sheet  used 
for  the  interrupted  crack  propagation  tests.  The  following  is  a  brief 
description  of  ‘'he  key  steps  in  a  typical  test  firing: 

a)  measure  the  main  ignitor  charge  and  shavings  and  assemble 
ignitor.  Measure  the  resistance  of  the  primer  and  short  the  primer 
circuit; 

b)  measure  and  record  the  crack  geometry  of  the  propellant  sample. 
Cover  all  surfaces,  except  crack  cavity,  with  flame  retardant  grease  to 
prevent  flame  spreading  over  those  surfaces.  Install  the  sample  and 
bolt  on  the  window  and  retainer; 

c)  install  pressure  transducers  and  calibrate  the  charge 
amplifiers.  Connect  amplifier  outputs  to  the  Physical  Data  recorder; 

d)  set  up  the  high-speed  camera  and  test  ingition  circuit; 

e)  put  charge  amplifiers  in  the  OPERATE  position,  reset  the 
Physical  Data  recorder  and  start  the  camera; 

f)  when  the  rupture  disk  bursts,  open  the  nitrogen  line; 

g)  turn  off  charge  amplifiers,  store  and  plot  the  pressure  traces 


and  disassemble  the  chamber  to  examine  the  recovered  sample. 


Chapter  3 


THEORETICAL  ANALYSIS 


3. 1  Crack  Sample  Geometry 

Figure  6  shows  the  damaged  sample  geometry  considered  in  this 
analysis.  From  this  geometry,  some  key  parameters  can  be  defined. 
First,  and  i  are  direct  geometric  factors  affecting  the  crack 
propagation  process,  and  the  number  of  cracks  n,  is  a  key  factor  in 
defining  the  degree  of  damage.11  The  affected  volume,  or  damaged-zone 
volume  is  defined  as 


V 

a 


"LDb 


(3) 


The  burning  surface  surface  area  in  the  affected  region  is  given  by 


Afa  =  2tr(n+1)LDb 


For  this  analysis,  a  specific  burning  surface  area  is  defined  as 


b  2  (n+1 ) 


s  V 


(4) 


(5) 


The  functional  forms  of  equations  (3)  and  (5)  are 

Va  =  Va(LD,  b)  (6) 

As  =  As(Ld,  n)  (7) 

3.2  Determination  of  Mechanical  Properties 

A  series  of  mechanical  tests  was  conducted  by  Harbert  and 
12 

Schapery  to  determine  the  mechanical  properties  of  the  two  propellants 


j*#  - -i  'l'i 
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(ARC  5051  and  ARC  4525)  under  non-burning  conditions. 


3.2.1  Relaxation  Modulus 

Uniaxial  constant  strain  tests  on  propellant  bar  specimens  were 
performed  at  various  temperatures.  Relaxation  moduli  versus  time  data 
are  plotted  for  each  temperature.  In  order  to  obtain  a  continuous  curve 
for  all  temperatures,  the  relaxation  moduli  are  plotted  versus  a  shifted 
time  where 

c  =  t/aT  (8) 

The  time  shift  factor  a^  was  determined  based  on  a  reference  temperature 
of  297  K  (i.e.,  a^.  =  1.0  at  297  K).  Figure  7  shows  the  functional 
relationship  between  temperature  and  time  shift  factor.  The  resulting 
master  relaxation  curve  for  each  propellant  is  shown  in  Figure  8. 


3.2.2  Crack  Propagation  Velocity 

Crack  growth  tests  were  conducted  on  strip  specimens  with  a  precut 
crack  (Fig.  9)  at  different  crosshead  rates.  The  crack  velocity  and 
strain  were  measured  directly  and  the  stress  intensity  was  calculated 
using 


K  =  o\J h/2(1-v  ) 


where  Poisson's  ratio  v=  0.5  and  a  is  the  stress  far  from  the  crack  tip 

exerted  by  the  crosshead. ^  Figure  10  shows  the  relationship  between 

crack  propagation  velocity  and  stress  intensity  factor  for  ARC  4525 

propellant.  It  can  be  seen  that  the  lines  of  constant  strain  are 

parallel  (i.e.,  independent  of  strain).  This  fact  suggests  the 

applicability  of  linear  viscoelastic  theory  to  characterize  the  crack 

14 

propagation  velocity.  The  results  of  this  theory  can  be  expressed  by 
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the  following  power  law 


— <K  )q 
aT  1 


In  this  equation  q  is  the  strain- independent  slope  obtained  from  Figure 
10,  and  C  is  a  material  parameter  which  varies  with  strain.  The 
functional  representation  of  equation  (10)  is  then 


Vp  =  Vp(Kr  e,  aT) 


3.2.3  Fracture  Energy  and  Related  Parameters 

From  an  idealized  crack  propagation  model  for  viscoelastic  media  as 
2  3-5 

used  by  Knauss  and  Schapery  ,  the  fracture  energy  (D  can  be  defined 
as  che  energy  required  to  elongate  a  segment  of  area  dA  to  the  point  of 
failure.  This  quantity  has  the  functional  relationship 


r  =  r(Kr  vp,  afm) 


where  is  the  maximum  stress  in  the  failure  zone.  The  stress 
intensity  factor  is  generally  a  function  of  the  geometry  of  the 
continuum  in  which  the  failure  zone  is  embedded  and  the  stress  in  that 
zone,  represented  by  chamber  pressure  in  this  application.  Thus 


KI  =  W  V  V  p) 


KT  =  K  (14) 

C 
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for  crack  initiation.  From  Langlois  and  Gonard,  the  extended  crack 
length  can  be  related  to  stress  intensity  factor  and  the  maximum  stress 
in  the  failure  zone 


j 


An  energy  balance  is  applied  to  the  control  volume  (Figure  11)  in 


order  to  relate  key  parameters  in  the  combustion  process.  The  basic 
assumptions  are: 

a)  The  heat  loss  from  the  control  volume  to  the  surroundings  is 
negligible  since  the  damage  occurs  in  an  extremely  short  time  interval 
(of  the  order  of  10  ms). 

b)  The  combustion  gases  of  the  driving  motor  are  considered  to 
follow  the  perfect  gas  law.  This  assumption  is  valid  since  the  gas 
temperature  is  high  (2000  K)  and  the  maximum  chamber  pressure  is  less 
than  35  MPa. 

c)  mQut  =  0  for  the  time  of  interest  since  the  initial  damage 
occurs  before  the  rupture  disk  bursts. 

d)  The  gas  phase  control  volume  can  be  treated  as  a  constant  since 
the  volume  increase  due  to  flow  work  and  surface  burning  is  extremely 
small  compared  to  the  total  control  volume  during  the  time  of  interest. 
The  conservation  of  energy  for  the  gas  phase  can  be  written  as 
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(4  F  •  C 

—  =  ft.  h.  .  +  ft  h  -  ft  ,  c  T  -  Q  -  P— rr- 
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(16) 


Using  the  basic  assumptions,  equation  (16)  reduces  to 


c  3P  .  .  .  . 

— 7  —  =  ft.  h.  .  +  m  h 

Y-1  3t  ign  t ,ign  p  p 


After  replacing  according  to 


(17) 


where  the  total  burning  surface  area  can  further  be  related  to  the 
crack  geometry 


Abt  =  Abt(Va’  As’  LD’  ao* 


Equation  (17)  then  takes  on  the  functional  form 


— -  =  — (V  ,  ft.  ,  h„  .  ,  r.,  1  ,  !  ,  L,,  p  ,  (1,  1  ,  h  )  (20) 

at  at  c  lgn  l-lg"  b  a  0  D  p  1  3  p 

The  burning  rate  under  a  pressure  transient  has  been  characterized 

.  „  .  16 
by  Kner  as 

<j>na  . 

r  (t)  =  r  [1  +  - B-r  (21  ) 

b  °  P(t)r2dt 

o 

where  rQ  is  the  steady  state  burning  rate  according  to  the  Saint-Robert 
17 

law,  ap  is  the  solid  phase  thermal  diffusivity,  and  n  and  $  are 

numerical  constants.  Since  the  thermal  diffusivities  were  the  same  for 
12 

both  propellants  it  also  was  treated  as  a  numerical  constant.  The 

instantaneous  burning  rate  can  therefore  be  related  as 

r.  =  r.  (P,  3P/3t)  (22) 

b  b 

Since  the  flow  from  the  driving  motor  is  choked  and  the  ignitor 
charge  is  made  of  the  same  propellant  as  the  sample,  the  propellant 
enthalpy  can  be  related  to  the  stagnation  enthalpy  of  the  igniter  gases 


hp  ”  Y-1ht,ign 


The  propellants  are  essentially  incompressible  and  of  equal 


density,  therefore 


< 


constant 


(24) 


From  evidence  obtained  in  crack  propagation  experiments,^  crack 
propagation  velocity  in  the  burning  solid  propellants  is  related  to  the 
following  parameters: 


V  =  V  (3p/3t, 
P  P 


r’  V 


aT> 


(25) 


3.4  Dimensional  Analysis 

Equations  (6),  (7),  (11 )— (15),  (20)  and  (22)-(25)  relate  the 
twenty-three  parameters  (y  was  assumed  to  be  a  numerical  constant  for 
the  combustion  product  gases)  affecting  the  crack  propagation  process. 
These  twelve  equations  may  be  reduced  to  one  equation  involving  twelve 
of  the  more  significant  parameters,  namely 

Vp  =  yavat.  v  .o,  y  k  y  *s.  y  y  y  aT>  <26> 


3.4.1  Arrangement  into  Meaningful  Dimensionless  Groups 

1  8 

The  Buckingham  Pi  Theorem  indicates  that  there  exist  nine 

independent  dimensionless  parameters  involved  in  the  crack  propagation 

process.  Three  of  the  key  parameters  are  already  dimensionless,  thus 

4^,  4^  and  a^  were  chosen  as  dimensionless  parameters  describing  the 

geometry  of  the  continuum  around  the  damaged  zone  and  the  effect  of 

initial  propellant  temperature  on  mechanical  properties.  In  order  to 

develop  additional  meaningful  groups,  the  parameters  which  describe  the 

degree  of  damage  (i.e.,  Ln,  V  ,  A  and  V  )  were  examined  individually. 

u  a  s  p 

A  desireable  dimensionless  parameter  defining  the  extended  crack 
length  was  the  percentage  increase  in  overall  crack  length.  This 
dimensionless  parameter  was  defined  as 


vs 
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^D  _  Length  of  extended  macrocrack 
aQ  ~  Length  of  the  initial  crack 


(27) 


From  the  parameters  V  and  A  ,  two  meaningful  dimensionless  groups 

3  S 

could  be  formed, 

V 

a  s  _  Rate  of  increase  of  crack  surface  area  (28) 

B  ”  2Lp  "  Rate  of  increases  of  cross-sectional 
area  of  the  affected  region 

ttc  =  T^Lp  -  1  =  n,  number  of  macrocracks  of  length  LD  (29) 

The  crack  propagation  velocity  was  normalized  to  the  normal 
regression  rate  of  the  surface  due  to  burning, 


■nn  = 


Crack  propagation  Velocity _ 

Burning  rate  of  the  main  propellant 


(31  ) 


The  remaining  dimensionless  parameters  relate  the  input  energy 


from  the  propellant  and  the  presurization  rate  to  the  critical  values 
for  crack  propagation  and  were  grouped  as  follows: 


TT 


E 


hppprbAsTrLd 
V  [K^  (1~v2)/E 

P  xc 


ref 
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Rate  of  thermal  energy  input 
from  propellant  product  gases  (32 ) 
Strain  energy  release  rate 
required  for  crack  propagation 


9P  ^  3/2  Rate  of  change  of  pressure  force  exerted 

_  dt  d _ on  the  crack  surface _  (33) 

F  Kj.  rfa  ”  (critical  stress  for  fracture)x(rate  of 

C  increase  of  flow  area  in  the  damaged  zone) 


The  additional  constants  such  as  density,  sample  thickness  and 
Poisson's  ratio  were  added  to  give  the  groups  their  respective  physical 
meaning . 
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3.^.2  Functional  Relationships  Between  Dimensionless  Parameters 
The  dimensionless  parameters  tt a,  tt^ ,  tt  ttd  characterize  the 
degree  of  damage  in  the  propellant  sample.  The  driving  force  causing 
this  damage  (i.e.,  the  initial  pressurization  rate  and  thermal  energy 
input  from  the  propellant  gases)  are  characterized  by  the  dimensionless 
parameters  and  tt^,.  These  parameters  were  used  together  with  ^ ^ 
and  a^,  for  constructing  the  following  emperical  correlations  to  predict 


the  crack  propagation  and  branching  processes: 


f  y  • 
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Chapter  4 
RESULTS 

4.  1  Determination  of  Dimensionless  Parameters 

Several  tests  were  conducted  using  the  aparatus  previously 
described.  From  the  pressure  data,  film  data  and  known  mechanical  and 
thermodynamic  properties,  values  for  each  of  the  parameters  could  be 
deduced . 

4.1.1  Pressure  Loading  in  the  Damaged  Zone 

Figure  12  shows  the  pressure-time  trace  at  three  key  locations  in 
a  typical  test.  Location  1  corresponds  to  the  chamber  free-volume  and 
location  2  corresponds  to  the  initial  crack  tip  location  and  the  3rd  is 
further  along  the  crack  propagation  direction.  The  first  curve  is  the 
driving  pressurization  rate  used  in  calculating  the  dimensionless 
parameters  since  it  is  generally  independent  of  the  initial  crack 
configuration.  The  pressure  trace  for  the  second  location  illustrates 
the  increased  pressurization  rate  and  peak  pressure  experienced  at  the 
tip  of  an  opening  crack.  This  may  be  caused  by  the  delayed  opening  of 
the  crack  and/or  the  increased  burning  rate  in  the  failure  zone.  As 
seen,  the  peak  pressure  in  this  case  was  13.8  MPa  (2000  psi)  when  the 
rupture  disk  burst.  The  chamber  pressure  then  dropped  rapidly,  and  the 
sample  was  recovered.  Figure  13  is  a  close-up  photograph  of  the 
damaged  zone  showing  the  bracnching  which  occurred. 

Figure  14  is  a  similar  pressure  trace  for  another  test.  The 
rupture  disk  burst  earlier  in  this  case  (7.0  MPa),  which  may  have 
contributed  to  the  reignition  and  burn-out  of  the  sample.  The  early 
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rupture  appeared  to  affect  the  peak  pressure  at  the  crack  tip  as  it 
occurred  during  the  depressurization  of  the  free  volume. 


4.1.2  Interpretation  of  the  Film  Data 

The  degree  of  damage  to  the  propellant  sample  was  characterized  by 
examining  the  film  record  of  the  event.  A  typical  film  record  of  a 
Mode  B,  single  crack  propagation  is  shown  in  Figure  15.  Table  I 
contains  a  frame-by-frame  decription  of  the  event.  The  instantaneous 
bright  flame  contours  caused  by  burning  in  the  crack  are  used  to  define 
the  crack  shape  and  length.  From  the  contour  (Figure  16),  the  crack 
tip  location  was  measured  as  a  function  of  time.  Thus  the  crack 
velocity,  extended  length  and  number  of  macrocracks  were  deduced.  In 
this  test  the  initial  crack  length  was  longer  than  in  previous  tests. 
Figure  17  shows  that  the  luminous  flame  front  moved  through  the 
initially  closed  region  and  the  crack  propagation  region  at  the  same 
rate.  The  sharp  reduction  in  propagation  rate  was  caused  by  the 
depressurization  of  the  chamber.  A  typical  instantaneous  bright  flame 
contour  for  a  branching  experiment  is  shown  in  Figure  18. 

4.1.3  Determining  Other  Key  Parameters 

In  order  to  calculate  each  dimensionless  parameter  the  values  of 

rb,  Kj  and  Eref  must  be  determined.  To  represent  the  normal 
C 

regression  of  the  propellant  surface  in  relation  to  the  crack 

propagation,  a  steady  state  burning  rate  was  calculated  using  Saint 
17 

Robert's  law 

rfe  =  aPn  (34) 

where  P  was  the  pressure  at  the  time  of  crack  initiation  and  a  and  n 
are  burning  rate  constants  given  in  Table  II. ^ 
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Table  I.  Film  Interpretation  of  Crack  Propagation  Process 


Picture  No.  Time  (ms) 


13-16 


17-28 


29-34 


0-0.6 


0.6-1. 2 


1.2-1. 6 


1 .6-2.8 


2. 8-3. 4 


Event 


Igniter  gases  fill  crack  cavity. 


Flame  penetrates  precut  crack. 


A  single  crack  propagates.  High 
pressure  gas  continues  to  flow 
into  crack  region. 


Crack  cavity  is  clearly  seen. 
Propagation  has  stopped. 


Rupture  disk  bursts  and  flame 
intensity  is  reduced. 


Propellant  sample  is  not  extinguished 
and  burns  completely. 
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Table  II.  Steady  State  Burning  Rate  Constants 
Type _ a  (mm/s/(atm)n  exponent  n 


ARC  4525 
ARC  5051 


0.5849 

0.8441 


0.5427 
0.561 1 


The  critical  stress  intensity  factor  K_  ,  was  calculated  using 

C 

equations  (9),  (14)  and  the  relationship 


o  =  Pcos^ 


(35) 


The  relaxation  modulus  was  taken  from  Figure  8  where  t  was  the 

time  of  crack  initiation. 


4.2  Relationships  Between  Dimensionless  Parameters 

Using  the  data  obtained  from  numerous  experiments,  several  sets  of 
values  for  the  nine  dimensionless  parameters  have  been  deduced  (Tabic 
III).  With  this  data,  the  functional  relationships  of  these  parameters 
were  obtained  statistically  using  a  multiple  regression  analysis 
program.  The  results  can  be  expressed  as 


,n„  0.107  0.181  0.321,  .  .  .2.49,  .  ,  .0.645 

irA  =  20.697  Tip.  aT  (sinijO  (simpo) 


^  ,n-T  0.065  0.254  -0.78,  .  ,  .-6.45,  .  .  .-0.985 

tt—  =  6.16x10  m.  tt_  a_  (simp.)  (smtp  ) 
d  fc.  r  1  l  o 


„  o0n  0.274  0.042  -0.021,  .  ,  .1.14,  .  .  .-0.045 

ttc  =  7.389  tte  Up  aT  (sini|j.)  (simpo) 


n  ,-4  0.217  0.076  0.478,  .  .  ,2.05,  .  ,-0.117 

iT—  =  10.9x10  rr-  tt_  a_  (simp.)  (smw  ) 

D  E  F  T  i  o 


(36) 


(37) 


(38) 


(39) 


It  is  important  to  note  that  although  the  exponents  of  and 
are  small  compared  to  the  other  exponents,  the  values  of  these  two 
parameters  are,  in  general,  large  compared  to  simp^,  sinipQ  and  a^,  and 
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therefore  have  a  dominant  effect. 


In  order  to  check  the  reliability  of  these  predictions,  equations 

(36)  through  (39)  are  plotted  with  the  values  obtained  from  test  data 

(Figures  19  through  22).  These  figures  show  that  the  predictions 

adequately  characterize  the  experimental  data.  The  multiple  regression 

2 

analysis  calculated  R  values  greater  than  89%  for  all  correlations. 

2 

R  is  defined  as  the  sample  multiple  coefficient  of  determination  given 
by 


R 


2 


~  2 

E(Y Y. ) 

1  1 

Z(Y . -  Y)2 
1 


(40) 


wnere  Y.  are  the  measured  values,  Y  is  the  mean  value  of  Y.  and  Y.  is 
1  11 

the  predicted  value  of  Y^. 


4. 3  Discussion  of  Results 

The  expressions  for  tt^,  tt^,  and  Up  characterize  the  damage  in 

the  sample  in  terms  of  the  extended  macrocrack  length,  in  rease  in 

burning  surface  area,  number  of  macrocracks,  and  crack  propagation 

velocity.  From  equations  (36)  through  (39),  the  effects  of  the 

operating  conditions  characterized  by  it  tt  ip  ,  1 p.  and  a„  on  the 

b  r  o  1  I 

sample  damage  can  be  analyzed  individually.  The  equations  show  that 
all  damage  parameters  increase  with  tt  and  tt  The  pressurization  rate 

L  r 

(tt p )  is  particularly  dominant  in  extending  the  crack  length  and 
increasing  the  specific  surface  area.  This  agrees  with  previous 
studies  which  showed  an  increase  in  damage  for  higher  rates  of 
thermal  energy  input  and  higher  pressurization  rates. 

Changes  in  p  change  the  confinement  of  the  damaged  zone  allowing 
varying  degrees  of  sample  deformation.  The  equations  predict  that  a 


greater  angle  will  increase  crack  growth  while  each  of  the  other 

parameters  will  decrease.  Under  these  conditions  the  angle  at  the 

crack  tip  is  likely  to  be  greater  resulting  in  a  higher  local  stress 

concentration.  Thus,  the  single  crack  propagation  mode  is  more 

probable.  With  a  longer,  single  crack,  equation  (37)  indicates  that 

these  two  effects  contribute  to  a  decrease  in  tt  Also,  higher  sample 

o 

deformation  work  implies  less  energy  is  available  to  propagate  the 
crack,  reducing  V  . 

The  effects  of  initial  propellant  temperature  on  the  crack 
propagation  process  can  be  very  important.  At  lower  temperatures 
(higher  values  of  a^)  the  propellant  becomes  more  brittle,  lowering  its 
fracture  energy.  Higher  crack  velocities  and  longer  cracks  are 
expected.  However,  the  range  of  initial  temperatures  in  this  study  was 
limited  by  the  test  conditions.  Broader  ranges  are  needed  to  more 
accurately  predict  its  effect,  i+k  was  essentially  unchanged  in  these 
tests  but  was  included  for  the  completeness  of  the  correlations. 

The  agreement  of  these  predictions  with  previous  observations  and 
basic  principles  in  fracture  mechanics  indicates  that  the  methodology 
for  their  development  is  valid.  Similar  correlations  for  different 
propellant  types  with  different  crack  configurations  could  be  developed 
using  this  method  to  determine  guidelines  for  the  safe  operation  of 
high  performance  rocket  motors. 
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Chapter  5 

SUMMARY  AND  CONCLUSIONS 


Two  types  of  composite  propellants  were  selected  for  this  study 
(ARC  5051  and  ARC  4525).  They  have  similar  composition  but  different 
oxidizer  particle  sizes  and  mechanical  properties.  Several  tests  were 
conducted  under  different  loading  conditions  and  sapmle  geometries. 

The  damage  to  the  sample  was  measured  in  terms  of  crack  length,  number 
of  cracks  and  crack  propagation  velocity.  Using  a  dimensional  analysis 
of  the  key  variables  involved  in  the  crack  propagation  phenomena,  a  set 
of  governing  dimensionless  parameters  were  established.  These 
parameters  can  provide  some  insight  into  the  cause  and  effect 
relationship  between  the  momentum  and  energy  input  to  a  propellant 
sample  with  existing  flaws  and  the  resultant  damage. 

Correlations  of  the  dimensionless  parameters  were  obtained  by 
performing  a  multiple  regression  analysis  on  the  test  data.  The 
driving  parameters,  tt^.  and  have  positive  relationships  with  the 
damage  parameters  Tig ,  ttc  and  1^.  This  indicates  that  higher 

momentum  and  energy  inputs  to  the  crack  cavity  will  cause  a  greater 
degree  of  damage.  The  geometry  and  initial  temperature  of  the 
propellant  sample  have  additional  effects  characterized  by  the 
correlations.  These  correlations  are  only  valid  for  the  small  ranges 
of  temperature  and  crack  cavity  angles  studied;  however,  additional 
tests  could  expand  the  applicability  and  validity  of  this  method. 

Enough  data  has  been  obtained  to  predict  the  strong  dependence  on  the 
rate  of  energy  input  and  pressurization  on  crack  propagation.  The  good 
agreement  with  test  data  and  reasonable  physical  interpretation 
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suggests  that  the  methodology  outlined  in  this  study  can  provide  a 
basis  for  future  development  of  similar  correlations  defining  safe 
operating  conditions  for  solid  propellant  rockets. 

Future  research  in  this  area  should  include: 

a)  development  of  correlations  for  high-energy  propellants  used  in 
actual  propulsion  systems  today,  including  high-elongation  propellants; 

b)  extension  of  the  present  analysis  to  include  propagation  and 
flame  spreading  along  a  motor  casing  due  to  an  initially  debonded  area, 
and; 

c)  integration  of  a  crack/debond  model  into  the  Rocket  Performance 
Prediction  Code. 
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APPENDIX  A 


TEST  PROCEDURES  AND  CHECKLIST  FOR  CRACK  PROPAGATION  TESTS 

A.  Igniter  Assembly: 

_  1.  Cut  igniter  charge  and  record  mass  and  propellant  type. 

_  2.  Select  a  primer  with  a  resistance  20  -  50  ohms  and  record. 

_  3.  Assemble  primer,  primer  holder  and  lexane  insulator. 

_  4.  Check  primer  resistance.  If  it  has  changed  significantly, 

check  leads  and  reassemble. 

_  5.  Place  desired  amount  of  propellant  shavings  in  and  around 

the  primer. 

_  6.  Assemble  remaining  igniter  pieces  with  the  yellow  marks 

aligned. 

_  7.  Check  igniter  circuit  continuity. 

_  8.  Short  igniter  leads. 

B.  Sample  preparation: 

_  1.  Measure  and  record  the  mass  and  dimensions  of  the  sample. 

__ _  2.  Paint  a  dot  on  the  front  surface  of  the  propellant  for 

focusing  purposes. 

_  3.  If  Spin  Physics  camera  will  be  used,  mark  the  crack  tip 

location  with  paint  and  measure  distance  between  crack  tip  and 
focusing  dot. 

C.  Chamber  Assembly: 

_  1.  Install  the  desired  inlet  lines  (except  on  igniter)  and 

check  for  unobstructed  flow. 

_  2.  Assure  sacrificial  and  main  windows  properly  fit  into  test 

chamber  and  clean  window  surfaces. 


_  3.  Spread  flame  retardant  grease  over  all  sample  surfaces 

(except  crack  surface)  and  install  in  test  chamber. 

_  4.  Insure  propellant  sample  will  not  extend  beyond  sample 

holders  when  compressed. 

_  5.  Install  main  0-ring  with  joined  ends  at  the  rear  of  the 

chamber.  Apply  permatex  to  O-ring  and  groove  . 

_  6.  Install  sacrificial  window  against  the  propellant  sample. 

_  7.  Install  the  main  window  into  the  chamber  window  cover. 

_  8.  Bolt  the  chamber  window  assembly  to  the  main  chamber, 

being  careful  to  evenly  compress  the  sample. 

_  9.  Measure  and  record  top  bursting  diaphragm  thickness  and 

install  diaphragm  and  hood  assembly  to  top  of  chamber.  Check  for 
proper  alignment  of  diaphragm  holder  to  allow  maximum  flow  area. 

-  10.  Check  side  rupture  disk  and  replace  if  necessary. 

_  11.  Install  igniter  assembly  to  bottom  of  chamber  being  sure 

to  align  the  yellow  marks.  Check  igniter  resistance  and  short  the 
igniter  circuit. 

_  12.  Connect  line  to  igniter.  If  not  used,  ensure  that 

igniter  ports  are  plugged. 

D.  Instrumentation: 

_  1.  Verify  that  all  ignition  circuits  are  in  the  off  position. 

_  2.  Install  pressure  transducers  and  record  their  respective 

locations,  sensitivities,  and  data  lines. 

_  3.  Connect  water  lines  and  check  for  leaks.  Turn  off  water 

supply. 

_  4.  Set  up  Hycam  and  install  dummy  film.  Insure  the  film 

speed  is  reduced  for  ignition  circuit  test. 
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5.  Connect  volt-meter  to  ignition  box  and  run  ignition  test. 


6.  Check  focus  of  the  Hycam. 


7.  Load  new  film,  depress  microswitch,  and  insure  triggering 


arm  operates  properly  as  the  camera  door  is  secured.  Reset  PPS 


dial . 


8.  Check  and  record  all  camera  settings. 


9.  Record  test  cell  temperature. 


10.  Setup  and  focus  Spin  Physics  camera  (if  used). 


11.  Remove  lights  used  for  focusing. 


12.  Install  exaust  fan  and  connect  to  Lights  No.  1  receptacle. 


13.  Check  data  line  connections  to  charge  amplifiers. 


14.  Set  input  voltages  and  calibrate  charge  amplifiers. 


15.  Connect  charge  amplifier  output  to  Physical  Data  recorder 


and  record  settings.  Turn  on  Physical  Data  unit  and  hit  reset 


button. 


E.  Test  firing  procedure: 


1.  Perform  final  check  of  all  connections. 


Hycam 


Spin  Physics 


Transducers 


Water  connections 


Exhaust  fan 


N2  lines 


Physical  Data 


Charge  Amps 


2.  Open  outside  door. 


3.  Secure  safety  key. 


4.  Connect  igniter  leads  to  igniter  box. 


5.  Check  igniter  resistance. 


6.  Close  and  secure  test  cell  door. 


7.  Turn  on  warning  light  and  sound  test  alarm. 


ft 
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_  8.  Verify  charge  amps  are  in  CHARGE  position. 

_  9.  Turn  on  igniter  power  and  charge  capacitor. 

_  10.  Turn  on  ignition  panel  switches  3  thru  6  and  event  switch. 

_  11.  Turn  off  safety  switch. 

_  12.  Unground  charge  amplifiers,  reset  Physical  Data  recorder. 

_  13.  Turn  on  Spin  Physics  Camera  and  wait  for  ready  light. 

_  14.  Turn  on  ignition  switch. 

_  15.  After  rupture  disk  bursts,  inject  N^. 

_  16.  When  propellant  has  stopped  burning,  stop  Spin  Physics 

camera. 

F.  Post  Firing  Procedures: 

_  2.  Ground  charge  amplifiers. 

_  3.  Turn  off  all  ignition  switches  and  turn  on  exhaust  fan 

switch. 

_  4.  Make  a  hard  copy  of  pressure  traces  on  X-Y  plotter. 

_  5.  Transfer  pressure  data  to  Nicolet  and  store  on  floppy 

disk. 

_  6.  Generate  calibration  signals  for  transducers  and  record  on 


X-Y  plotter . 


DATA  SHEET  FOR  CRACK  PROPAGATION  EXPERIMENTS 

! 

CRACK  PROPAGATION  TEST 

NO. 

) 

DATE 

\ 

\ 

i 

Participants:  ,  , 

/  f 

< 

TEMPERATURE:  T(rocm)  C 

TEcell)  C 

< 

* 

CRACK  GEOMETRY  AND  SAMPLE  SPEC.: 

IGNITER: 

i 

1 

Propellant  type 

Primer 

! 

► 

Configuration 

Propellant  type 

i 

Length 

Dimensions 

Max  gap  width 

Mass 

I 

i 

Inner  taper  angle 

Resistance 

1 

Taper  heigth 

I 

S 

Thickness 

FLYING  PIN  IGNITER: 

; 

Weight 

Primer  resistance 

Fuse  Wire  Location 

Propellant  weight 

i 

Propellant  type 

* 

CAMERA: 

Diaphragm  thickness 

\ 

i 

Type  of  Lens 

Diaphragm  diameter 

! 

W 

Framing  Rate  Frames/s 

Diaphragm  material 

« 

* 

Cptical  Head:  1/4  ,  1/2  ,  frame 

r* 

i 

PPS  Dial 

RUPTURE  DISK  ACTUATOR: 

", 

Multiplier 

Primer  resistance 

.• 

Time  Marker  Pulses/s 

Propellant  weight 

» 

Remote  firing  at  ft 

Propellant  type 

> 

t 

V 

Distance  of  focus  ft 

V 

l 

Brake  Dial 

FILM: 

» 

Servo 

ASA  No. 

> 

§ 

s. 

V 

H/L  speed  reducer 

Length 

F.  Stoo: 

Color  ,  B  &  W 

► 

r 

/ 

■r 
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APPENDIX  B 


ERROR  ANALYSIS 

In  the  multiple  regression  of  the  data  in  Table  II,  the 
exponential  relationships  assumed  for  the  dimensionless  parameters 
[equations  (36)  through  (39)]  were  rearranged  into  a  linear 
relationship  based  on  the  natural  logarithm  of  each  value.  Using  the 
basic  properties  of  logrithmic  representations,  the  uncertainty  in  each 
dimensionless  quantity  can  be  expressed  as 

ln(x.)  =  ln(x)  +  e(x)/x  (41) 

where  x.  is  the  actual  value,  x  is  the  measured  value  and  e(x)  is  the 

l 

20 

error  interval  associated  with  that  measurement.  The  method  of  least 
square  analysis  used  on  the  measured  data  will  give  more  accurate 
results  when 

ln(x)  »  e(x)/x  (42) 

An  expression  for  the  relative  error  in  each  dimensionless  parameter  is 

r^  =  e(x)/(xln(x) )  (43) 

Table  IV  contains  estimations  of  the  various  errors  involved  in 
measuring  each  of  the  variables  in  the  experimental  setup.  Equation 
(43)  can  then  be  evaluated  for  each  of  the  dimensionless  parameters. 
This  equation  indicates  that  for  values  of  x  which  approach  unity,  the 
relative  error  could  be  significant.  Since  the  error  of  (n)  is 
assumed  to  be  zero,  its  correlation  is  not  affected.  However,  does 

have  a  finite  error  associated  with  it.  For  the  worst  case  of 
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=  1.1,  =  0.095  or  9.5%  relative  error.  This  could  be  improved  by 

referencing  a^,  to  a  temperature  outside  the  normal  operating  range. 

Using  equation  (43)  and  the  values  listed  in  Table  4,  the  relative 
error  was  calculated  for  each  dimensionless  parameter  based  on  its 
average  value. 


7 r,  =  L^/a 
A  Do 


Assume  e(L^)  <<  e(L^)  therefore  e(7T^)/TT^  =  e(L^)/LD.  Based  on  the 

average  value,  tt,  =  0.326,  r  =  9.9%. 

A  TT. 

A 

Using  this  method,  the  following  values  of  rx  are  obtained: 


r  =  9.9 % 
r  =  3.9% 
r  =  0.0% 
r  =  0.3% 
r  =1.7? 
r  =  3.4% 


r  .  ,  =  5.9% 

sln^i 

r  =  8.2% 

CO  Sill 
To 

r  =  2.1% 

T 


vV-'v'S*'-'.' 


Table  IV.  Estimation  of  Measurement  Errors 


Physical 

Variable 


Source  and  TvDe  of  Error 


Error 

Interval  (  +  /-) 


Random  error  in  the  measurement  and 
calculation  of  the  inner  taper  angle. 

The  systematic  error  of  symmetry 


Random  error  in  the  measurement  of  the 
room  temperature 


0.5  K 


Systematic  error  in  the  difference  between 
room  temp,  and  initial  propellant  temp. 

Random  and  systematic  error  in  converting 
temperature  to  a_  (note:  in  the  temperature 
range  of  interest  da^./dT  is  near  zero) 

Accuracy  of  the  a^  data  (note:  aT  is 
defined  near  operating  temperatures) 

Since  n  is  a  discrete  quantity,  it  is  assumed 
to  have  no  random  error 


3.0  K 


Random  error  of  the  measuring  instrument 


0.25  mm 


Random  error  in  determining  the  bright  flame 
front  tip  location 


2.0  mm 


Systematic  error  of  correlation  between  the 
flame  front  and  actual  tip  location 
(note:  assumed  to  be  small) 


unknown 


Random  error  of  measuring  the  change  of 
flame  front  location  with  respect  to  time 
(due  to  the  averaging  effect) 


Systematic  error  of  the  installation,  wear, 
calibration  and  recording  of  pressure  data 


Since  dp/dt  is  calculated  using  a  very  large 
dp  and  t  is  very  accurate  the  error  should  not 
be  greater  than  in  P 


5% 


